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Sciences, 235B Animal Science/Vet. Med. BIdg., University of Minnesota, 1988 Fitch Avenue, St. Paul, MN 55108, USA. 

Toll-like receptors (TLRs), as innate immunity sensors, play critical roles in immune responses. Six SNPs of 
TLRS, TLR7, and TLRS were genotyped to determine their associations with systemic lupus erythematosus 
(SLE) and clinical manifestations of SLE. TLR7 SNP rs3853839 was independently associated with SLE 
susceptibility in females (G vs. C: p - 0.0051). TLR7 rs3853839-G (G vs. C: p - 0.0100) and TLRS 
rs3764880-G (recessive model: p — 0.0173; additive model: p — 0.0161) were associated with pericardial 
effusion in females relative to healthy females. Anti-SSA positive cases were more likely to have the 
dominant TLR7 rsl79010-T allele than normal controls (p - 0.0435). TLR3 rs3775296-T was associated 
with photosensitivity (p - 0.0020) and anemia (p - 0.0082). The "G-G" haplotype of TLR7 rs3853839 and 
TLR8 rs3764880 increased risk of SLE in females (age adjusted p — 0.0032). These findings suggest that TLR 
variations that modify gene expression affect risk for SLE susceptibility, clinical phenotype development, 
and production of autoantibodies. 

System lupus erythematosus (SLE) is a systemic autoimmune disease characterized by widespread loss of 
immune tolerance to self- antigens. Genetic and environmental factors contribute to the development of 
SLE, and patients typically experience alternating periods of flare-up and remission'. The generation of 
numerous autoantibodies that react with self-nuclear and -cytoplasmic antigens is associated with the dysfunc- 
tion of multiple organ systems^ ^\ The genetic transmission and patterns of inheritance of SLE have not yet been 
elucidated. In particular, it has been difficult to identify specific genetic polymorphisms associated with SLE due 
to the presence of polygenic inheritance, extensive genetic heterogeneity, the small size of most genetic studies, 
and the low disease prevalence''. Genome-wide association studies (GWASs) of various populations have 
identified several common immune response pathways and the presence of genetic variants in some ethnic 
groups that are involved in the pathogenesis of SLE* However, genetic dissection of SLE and other autoimmune 
diseases is difficult because these are complex diseases that involve alterations of multiple biologic pathways'*. 
The advent of modern genomics and the availability of new technologies have made it possible to fine-map 
candidate genes for SLE and other specific diseases based on knowledge of gene map position and functional 
relevance'^^'". Such studies may help to identify the roles of candidate genes for SLE and the roles of different 
genes in the expression of different clinical manifestations*"'^". 

Previous studies have indicated that innate pattern recognition receptors, such as Toll-like receptors (TLRs), 
play important roles in the development of autoimmunity. TLR proteins are localized on the cell surface or in 
endosomes, and play critical roles in innate immune responses against different pathogens""'"'^. Internalized 
nucleic acid immune complexes act as endogenous ligands that activate intracellular TLRs, and these initiate 
several signaling pathways that lead to increased production of type I interferons (IFNs) in plasmacytoid dendritic 
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Table 1 | Clinical characteristics of the SLE study population (n 


= 795) 
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= 795) 
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cells (pDCs)^^ Increased production of type I IFNs increases apop- 
tosis, neutrophil cell death via neutrophil extracellular trap 
(NETosis), innate immune signaling, and viral infection-induced 
autoimmunity^'"^'. Aberrant stimulation of the innate immune sys- 
tem through intracellular TLRs may lead to hyperactive immune 
responses and contribute to the pathogenesis of SLE"*" 

TLR3, TLR7, and TLRS are primarily associated with endosomal 
membranes and they recognize microbial nucleic acids. TLR3 binds 
to viral double-stranded RNA (dsRNA), and induces antiviral 
immune responses by promoting the production of type I IFN and 
pro-inflammatory cytokines. TLR7 binds to single-stranded RNA 
(ssRNA) from RNA viruses, and triggers pDCs to produce type I 
IFN. TLRS is phylogeneticaUy related to TLR7 and also mediates 
recognition of viral ssRNA. Thus, these 3 TLRs are responsible for 
pathogen clearance, antigen recognition, and induction of cytokine 
production'°•''■'^ 

Genetic and hormonal factors appear to partially explain the 
female predominance of SLE'^. In particular, TLR7 and TLRS are 
on the X chromosome (Xp22.2) and play critical roles in innate 
immunity and inflammatory responses^". Studies of animal models 
indicated that TLR7 gene dosage significantly affected the hyperac- 
tivity of B cells'"' ■". Mice with a Y-linked autoimmune accelerator 
(Yaa) that carries an extra copy of TLR7 develop autoimmunity to 
RNA-associated autoantigens, but mice with a low TLR7 copy num- 
ber require additional susceptibility loci to develop autoimmu- 
nity^'"'". In humans, males with Klinefelter syndrome carry an 
extra X chromosome (47, XXY) and are more likely to develop 
SLE; females with Turner syndrome lack one X chromosome (45, 
X) and are less likely to develop SLE''^ ''\ In addition, a study of a 
Mexican population indicated that increased TLR7 copy number 
correlated with TLR7 mRNA levels and susceptibility to child- 
hood-onset SLE'''' although TLR7 copy number variations (CNVs) 
are infrequent in human SLE''^. Other studies indicated that 
increased expression of TLRs in peripheral blood mononuclear cells 
and lymphocytes led increased IFN-cx expression in SLE patients'"" 

Previous research indicated that a functional TLR7 SNP 
(rs3S53839-G > C) affects TLR7 expression by modulation of 
microRNA-314S (miR-314S)'", and that two other intron SNPs 
(rs5935436 and rsl79010) were associated with increased SLE 



susceptibility''^''". The TLRS SNP rs37648S0 is a functional poly- 
morphism that affects TLRS transcription and translations of 
TLRS isoforms, which leads to the activation difference of NF-kB, 
is also associated with SLE susceptibility''^''""^^. A study of cell cul- 
tures indicated that TLR3 rs3775296 (in the promoter) and 
rs3775291 (in exon 4) affected TLR3 cell surface expression and 
localization, and subsequently influenced NF-kB cascade induction 
although the effect of these mutations has not been assessed in clin- 
ical studies^^. The present study examined the associations of the 
above mentioned 6 SNPs (3 in TLR7, 1 in TLRS, 2 in TLR3) with 
SLE and with specific clinical manifestations of SLE. 

Results 

We investigated the role of 6 SNPs in the susceptibility to SLE 
(rs3775296 and rs3775291 from TLR3; rs5935436, rsl79010, and 
rs3S53S39 from TLR7; and rs37648S0 from TLRS) by examination 
of 795 SLE patients (6S males and 727 females) and 1162 healthy 
controls (513 males and 649 females). The average age of cases was 
30.71 years (SD = 11.62, range: S to 77 years), 8.55% were males 
(31.46 ± 12.52 years-old), and 91.45% were females (30.64 ± 11.56 
years-old). The healthy controls had a mean age of 40.24 years (SD = 
lO.SS, range lS-64 years), with similar average ages of males and 
females (40.26 ± 9.26 years vs. 40.23 ± 12.02 years). The age differ- 
ence of SLE patients and healthy controls was statistically significant 
for both males and females, so we adjusted for age in the subsequent 
association analysis. Table 1 shows the clinical characteristics of the 
795 SLE patients, with independent statistical analyses for males and 
females. In the female cases and controls, there were no deviations 
from HWE in the six candidate SNPs. 

Table 2 shows the single-locus associations between the six can- 
didate SNPs and susceptibility to SLE. In males, there were no sig- 
nificant case-control associations after the false discovery rate (FDR) 
correction, possibly due to the small sample size. In females, there 
was a significant allelic association between SLE and rs3S53S39 in 
TLR7 (PpDR = 0.013, OR = 1.3S, 95% CI = 1.12-1.69). Moreover, 
the risk allele G had a recessive effect in females (GG vs. GC + CC: 
PpDR = 0.017, OR = 1.44, 95% CI = 1.14-1.83), suggesting that 
rs3S53S39 plays a role in development of SLE. We also evaluated 
the independent contributions of the 6 candidate SNPs (adjusted for 
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Table 3 | Associations of TLR SNPs with SLE cl inical phenotypes and autoantibodies in male SLE patients 
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"The former allele is the risk allele. 
''RAF: Risk allele frequency. 

"The p-value was calculated from chi-square tests or Fisher's exact tests when appropriate. 
''Pfdr'- The p-value was adjusted after the FDR correction. 



age) to SLE risk in females, and performed a multivariate logistic 
regression analysis that included the significant SNPs with the same 
genetic models. The rs3853839 SNP was the only significant suscept- 
ibility marker among the 6 examined SNPs (G V5. C:p = 0.0051, OR 
= 1.42,95% CI = 1.11-1.82). 

TLR SNP polymorphisms affected SLE phenotype and production 
of autoantibodies. Patients with SLE are present with heterogeneous 
clinical features, and have significant variations in the severity, 
nature, and spectrum of clinical involvement. Thus, we examined 
the effect of TLR SNP polymorphisms on SLE clinical parameters 
and phenotypes. Based on the clinical characteristics of males and 
females (Table 1), we initially performed two comparisons: (i) allele 
frequencies of SLE patients with each characteristic (" + " in Tables 
S1-S3) and SLE patients without the characteristic (" — " in Tables 
S1-S3); and (ii) allele frequencies of SLE patients with each 
characteristic and the normal controls ("normal" in Tables S1-S3). 
In males, there were no significant associations of the 6 candidate 
SNPs with SLE clinical manifestations after the FDR correction 
(Table 3, with more details in Tables SI and S3). Again, this may 
result from the small sample size. 

In females, TLR7 rs3853839 G risk allele was associated with sev- 
eral clinical manifestations of SLE (Table 4). The comparison of 
phenotype-positive cases and normal controls indicated significant 
associations of this allele with oral ulcer, arthritis, malar rash, photo- 
sensitivity, pericardial effusion, depressed complement, anti- 
dsDNA, anti-Sm, and anti-SSA (PpoR < 0.05 for all comparisons). 
Multivariate logistic regression analysis of the independent asso- 
ciation of each clinical characteristic with rs3853839 indicated that 
only pericardial effusion was significantly associated with the TLR7 
rs3853839 G risk allele. In particular, cases with pericardial effusion 
were more likely to have the G risk allele than normal controls {p = 
0.0100, OR = 2.82, 95% CI = 1.28-6.19). Table S2 shows that anti- 
SSA was associated with the rs 1790 10 T risk allele (T vs. C: Pfdr = 
0.0366, OR = 1.47, 95% CI = 1.11-1.94; CT + TT vs. CO. Pfdr = 
0.0037, OR = 1.93, 95% CI = 1.33-2.81). Multivariate logistic regres- 
sion analysis also showed that cases who were anti-SSA positive were 
more likely to carry the dominant rs 1790 10 T risk allele than normal 
controls (p = 0.0435, OR = 1.33, 95% CI = 1.01-1.75), but cases who 
were anti-SSA negative had the opposite tendency {p = 0.0299, OR = 
0.69, 95% CI = 0.49-0.96). 



Although rs3764880 in TLR 8 was not associated with SLE sus- 
ceptibility, the risk allele G was associated with oral ulcer with sig- 
nificant additive effects (Pfdr = 0.0232: OR = 1.66, 95% CI = 1.18- 
2.35, G vs. A: Pfdr = 0.0300, OR = 1.65, 95% CI = 1.17-2.32). 
However, incorporation of clinical variables with nominal p values 
below 0.05 into the multivariate logistic regression model indicated 
that pericardial effusion was the only clinical characteristic assoc- 
iated with rs3764880. In particular, pericardial effusion-positive 
cases were more likely to carry the rs3764880 G risk allele than 
normal controls, either with a recessive effect (p = 0.0173, OR = 
2.91, 95% CI = 1.21-7.02) or with an additive effect (p = 0.0161, OR 
= 2.96, 95% CI = 1.22-7.14). 

Analysis of the association of TLR3 SNPs with SLE clinical char- 
acteristics indicated no significant associations for males (Tables 3 
and S3). In females, comparison of phenotype-positive cases with 
phenotype-negative cases and of phenotype-positive cases with nor- 
mal controls indicated positive associations of rs3775296 with ane- 
mia (Table 4, TT vs. GG + GT: Pfdr = 0.0244, OR = 2.41, 95% CI = 
1.33-4.39; TT vs. GG + GT: Pfdr = 0.0373, OR = 2.11, 95% CI = 
1.23-3.65, respectively). Multivariate logistic regression analysis 
indicated that photosensitivity-negative cases and anemia- negative 
cases were associated with recessive rs3775296 T risk allele relative to 
normal controls (photosensitivity: p = 0.0020, OR = 2.40, 95% CI = 
1.38-4.19; anemia: p = 0.0082, OR = 0.44, 95% CI = 0.24-0.81). 

TLR haplotypes were associated with SLE susceptibility in females. 

In Figure 1, pair- wise LD measures r^ and D' for the four TLR7 and 
TLR8 SNPs rs3853839, rs5935436, rsl79010 and rs3764880 in the 
female healthy controls were presented. We estimated the frequ- 
encies of haplotypes formed by these 4 SNPs in females to 
determine the associations of different haplotypes with suscepti- 
bility to SLE. Table 5 shows the 4 inferred haplotypes that had 
estimated frequencies more than 5%. Performance of haplotype- 
trait association tests in female SLE patients and healthy controls 
indicated that the "C-C-C-G" haplotype protected against develop- 
ment of SLE (5.78% in cases vs. 8.55% in controls; permutation p- 
value = 0.0033, OR = 0.60, 95% CI = 0.44-0.82 after age 
adjustment). In Table 6 shows the estimates for haplotypes from 
rs3853839 and rs3764880 SNPs only. These results show that the 
effect of "C-G" on SLE susceptibility remained significant when we 
controlled for age (7.08% in cases vs. 8.39% in controls; permutation 
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Figure 1 | Pair- wise linkage disequilibrium patterns with r^ (left side) and D' (right side) measures of the four SNPs in TLR7 and TLR8 at Xp22.3 for 
649 healthy female controls. 



p-value = 0.0073, OR = 0.62, 95% CI = 0.46-0.83). Nevertheless, 
females with "G-G" haplotypes of rs3853839 and rs3764880 were 
significantly more likely to develop SLE (age adjusted p-value = 
0.0032, OR = 1.32, 95% CI = 1.10-1.60). 

Discussion 

Previous research with large samples indicated that the functional 
TLR7 SNP rs3853839-G > C was significantly associated with SLE 
in East Asians, especially in males*^. The current study identified the 
rs3853839-G allele of TLR7 as the main susceptibility marker in 
female SLE patients from Taiwan based on single-locus and multi- 
variate logistic regression analyses. This confirms the role of this SNP 
in the pathogenesis of SLE. The genetic effect of TLR7 on SLE appears 
to vary among different ethnic groups''^'^"'""'''. In particular, a previous 
study indicated that the TLR7 intronic SNPs rsl79019 and rsl79010 
were associated with SLE, independent of the 3' UTR SNP rs3853839 
in Japanese females^''. Other research indicated that the TLR7 SNP 
rsl79008 was not associated with SLE in a European population, but 
was significantly associated with SLE in BrazUians^^'^''. Recently, Deng 
et al. conducted a large trans-ancestral fine-mapping of European 
Americans, Afi^ican Americans, and Amerindian/Hispanics and iden- 
tified rs3853839-G as the only genetic risk variant for SLE in the 
TLR7-TLR8 region, although they did not confirm the male specific 
association. Notably, rs3853839-G appears to increase risk for SLE in 
different populations, although there are different frequencies in dif- 
ferent ethnic groups. This highlights the critical role of elevated TLR7 
expression in the pathogenesis of SLE, which in this case is mediated 
by slower mRNA degradation due to miR-3148 expression""". 



TLR7 and TLR8 contribute to antigen recognition and antibody 
production in the pathogenesis of SLE. TLR ligands, which are pre- 
sent in viruses and virus-like particles (VLPs), can directly stimulate 
B cells, and B cell responses are integrated with dual antigen-specific 
B cell receptor (BCR) and TLR engagement^^'^". B cells with up- 
regulated MYD88 expression become more responsive to TLR 
ligands and promote class -switch recombination^'. The inactivation 
or overexpression of genes that encode TLRs or of molecules that 
alter TLR signaling provides a bridge between the innate and adapt- 
ive immune systems, and this is critical to the presence of B cells 
defects in the pathogenesis of SLE'""''''. Moreover, the recognition of 
endogenous RNA-containing antigens by TLR7/3 may trigger auto- 
reactive B cells in the germinal center, and this is accompanied by 
the suppression of T regulatory cells, leading to disruption of self- 
tolerance*^ In this regard, numerous innate and adaptive immun- 
ity related genes involving IFN-alpha mediated signature pathway as 
well as T and B cells activation signaling pathways participate in the 
SLE pathogenesis''\ The present study found that TLR7 rs3853839 G 
risk allele was associated with several clinical manifestations of SLE, 
including oral ulcer, arthritis, malar rash, photosensitivity, pericar- 
dial effusion, depressed complement, and anti-dsDNA, anti-Sm, and 
anti-SSA autoantibodies. The TLR7 rsl79010 T risk allele was also 
associated with anti-SSA autoantibodies. These findings suggest that 
TLR7 may play a key role in autoantibody production because it 
increases B-cell sensitivity to RNA-containing autoantigens in the 
development of systemic autoimmunity. However, given the limited 
sample sizes of our stratified phenotype groups, this finding requires 
replication by larger future studies. 



Table 5 | Haplotype analysis of TLR7 and TLRS SNPs in female 


SLE patients 


and healthy controls 






Haplolype (rs3853839-rs5935436- 
rsl79010-rs3764880) 


Estimated Frequency {%) 


Permutation 


Logistic regression 


Logistic reg 


ression adjusted by age 


All SLE Normal 


p* 


P OR (95% CI) 


P 


OR (95% CI) 


G-C-C-G 
G-C-T-G 
C-C-C-A 
C-C-C-G 


44.00 45.32 42.53 
27.47 28.27 26.58 
10.54 9.56 11.63 
7.09 5.78 8.55 


0.1571 
0.2137 
0.0742 
0.0033 


0.1405 1.12(0.96-1.30) 
0.3229 1.09(0.92-1.29) 
0.0768 0.80(0.63-1.02) 
0.0049 0.66 (0.49-0.88) 


0.0849 
0.6663 
0.1279 
0.0016 


1.15 (0.98-1.36) 
1.04 (0.87-1.25) 
0.82 (0.63-1 .06) 
0.60 (0.44-0.82) 


*p-values for estimated haplotypes were generate 


d from 1 0,000 permutations using the expectation-maximization (EM) algorithm. 
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Table 6 | Haplotype analysis of the TLR7 SNP rs3853839 and the TLR8 SNP rs3764880 in female SLE patients and healthy controls 

Estimated Frequency {%) Permutation Logistic regression Logistic regression adjusted by age 



Estimated Haplotype 
(rs3853839-rs3764880) 



All 



SLE 



Normal 



p* 



OR (95% CI) 



OR (95% CI) 



G-G 
C-A 
C-G 
G-A 



75.00 
1 1.01 
8.39 
5.60 



77.37 
10.25 
7.08 
5.30 



72.34 
1 1.86 
9.86 
5.93 



0.0015 
0.2837 
0.0073 
0.3297 



0.0024 
0.1756 
0.0090 
0.4677 



1.31 (1.10-1.55) 
0.85 (0.67-1.08) 
0.70 (0.53-0.91) 
0.89 (0.64-1.23) 



0.0032 
0.2996 
0.0014 
0.7796 



1.32 (1.10-1.60) 
0.87 (0.68-1.13) 
0.62 (0.46-0.83) 
0.95 (0.67-1 .35) 



*p-values for estimated haplotypes were generated from 1 0,000 permutations using the EM algorithn" 



We also observed that the non-synonymous TLR8 SNP 
rs3764880-G allele was a risk factor for oral ulcer and pericardial 
effusion, with significant additive effects. Interestingly, previous 
research indicated that the TLR8 rs3764880-G allele protected 
against tissue damage in active tuberculosis and predicted a slower 
disease course in patients with HIV infections^^'^^. The TLR8 SNP 
rs3764880 alters the ATG start codon of TLR8 isoform B into a GTG. 
A methionine at position 4 of isoform B is used as the start codon for 
the TL_RS-rs3764880G allele, resulting in a truncated TLR8 isoform B 
with a shorter signal peptide ( 1 038 residues for the TLR8- 
rs3764880G allele vs. 1041 residues for the TLJ?S-rs3764880A allele). 
Immune cells carrying the TL_RS-rs3764880G allele had augmented 
TNFof-responses, but decreased translation of truncated TLR8 iso- 
form B and NF-kB production relative to those carrying the TLR8- 
rs3764880A allele'"■^^ Therefore, TLR8 appears to have important 
roles in autoimmune diseases and in response to infections. 

Previous research indicated that TLR7 and TLR8 have closely 
related functions in immune responses. We observed that the TLR7 
and TLR8 SNP haplotype rs5935436-C/rsl79010-C/rs3853839-C/ 
rs3764880-G protected against development of SLE, but the effect of 
rs3853839-C/rs3764880-G on SLE susceptibility remained significant 
when controlled for age. On the other hand, females with G-G hap- 
lotypes of rs3853839 and rs3764880 were significantly more likely to 
develop SLE. These results indicate that TLR8 may play a regulatory 
role in TLR7 function in innate immunity, similar to that documented 
in mice^^. 

TLR3 recognizes dsRNA, and its elevated level in human SLE 
peripheral blood cells suggests that it may play a key role in IFN 
signature gene activation^^''*^"**^. In addition, previous research indi- 
cated that TLR3 aggravated lupus nephritis in lupus-prone mice^**. 
We observed that TLR3 rs3775296-T risk allele had recessive effects 
on photosensitivity and anemia-negative SLE patients. Thus, our 
data indicate that TLR3 may play a role in the development of dif- 
ferent SLE phenotypes, and in antiviral responses that trigger 
expression of pro-inflammatory genes. 

In conclusion, the present study confirmed that certain functional 
TLR7 and TLR8 variations, particularly TLR7 rs3853839-G > C, 
modify gene expression and increase the risk for SLE, development 
of certain SLE phenotypes, and production of autoantibodies. These 
results suggest that TLR7/8 genetic variations are potential biomar- 
kers for prediction of SLE phenotypes, and also have implications for 
the development of therapeutic measures that may prevent various 
pathological conditions that are characteristic of SLE. 

Methods 

Characteristics of the study populations. All patients were recruited from the clinics 
of Chang Gang Memorial Hospital, and rheumatology specialists confirmed that all 
patients fulfilled the 1982 and 1997 American College of Rheumatology (ACR) 
diagnostic criteria for SLE^^. For the purpose of this study, healthy controls were 
selected following a questionnaire to assure that they did not have any autoimmune 
disease phenotypes. This study was approved by the ethics committees of Chang 
Gung Memorial Hospital, and all patients provided written informed consent 
according to the Declaration of Helsinki. 

Genomic DNA extraction. Genomic DNA was extracted from EDTA-anticoagulated 
peripheral blood using the Purgene DNA isolation kit as described previously™. 



SNPs genotyping assays. Validated Applied Biosystem TaqMan SNP assays were 
used for genotype determination of TLRs according the vendor's instructions (Life 
Technologies, Grand Island, NY). The TaqMan allele discrimination assays were 
performed on an ABI ViiA 7 Real-time System (Life Technology) and probes were 
labeled with a fluorescent dye (FAM and VIC). 

Statistical analysis. The functional candidate gene approach was used to perform a 
case-control association study by examination of 795 SLE patients (68 males and 727 
females) and 1162 healthy controls (513 males and 649 females). The Hardy- 
Weinberg equilibrium (HWE) was examined for the 6 selected SNPs using chi-square 
tests in females with and without SLE. Differences of allele frequencies in each of the 6 
SNPs were separately assessed in males and females to investigate the single-locus 
associations. Additionally, for females, the significance of differences in genotype 
frequencies were also evaluated, and dominant and recessive models were tested for 
each SNP. Thep-values, odds ratios (ORs), and 95% confidence intervals (CIs) were 
then calculated based on the risk allele identified. TLR 8 and TLR7 are located on the 
X chromosome (Xp22.2), so female and male data were analyzed separately. Meta- 
analysis with generation of meta p-values using the Cochran-Mantel-Haenzsel 
(CMH) method was used to assess the allelic associations of the 6 SNPs in all samples. 

We initially stratified clinical phenotypes according to diagnostic criteria to 
investigate the association of each SNP with SLE clinical manifestations. SLE patients 
with the phenotype under investigation were classed as cases, those without this 
phenotype as " — " cases, and healthy controls as "normal". Then, allele frequencies 
were compared for cases and " — " cases, and for cases and "normal" con- 
trols. In females, additional analyses also considered additive, dominant, and 
recessive effects for the risk allele of each SNP in order to assess genotype- phenotype 
associations. To investigate the independent association of SLE clinical characteristics 
with the six SNPs, multivariate logistic regression analysis was subsequently carried 
out to identify the independent statistical association of significant phenotypes 
(identified above) with the six SNPs. The additive, dominant and recessive allele 
effects for each SNP were modeled as the response variable and the three categories of 
individuals: "-F" cases, " — " cases and "Normal" pertaining to each clinical char- 
acteristic were used as the independent variables. 

Linkage disequilibrium (LD) patterns of neighboring SNPs on the same chro- 
mosome were analyzed by Haploview 4.2 (Broad Institute, Cambridge, MA, USA; 
http://www.broad.mit.edu/mpg/haploview). Haplotype information was inferred 
and frequencies were estimated using the HAPLOTYPE procedure in SAS 9.2 (SAS 
Institute, Cary, NC). Differences in haplotype frequencies were assessed for SLE cases 
and controls were separately assessed in males and females. The permutation (N — 
10,000) p-values of each haplotype were calculated using the expectation- max- 
imization (EM) algorithm, conditional on the other haplotypes, to evaluate the 
independent association of each category of haplotypes. Benjamini and Hochberg's 
linear step-up method in the SAS MULTTEST procedure was used to account for 
multiple testing^\ The False Discovery Rate (FDR)-adjustedp-values are defined in a 
step-up fashion, with less conservative multipliers and control. A corrected p-value 
(Pfdr) less than 0.05 indicated statistical significance. 
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